A note on Coulomb scattering amplitude 



m 
o 
o 



X 



Zafar Ahmed 

Nuclear Physics Division, Bhabha Atomic Research Centre, Bombay 400 085 
e-mail:zahmed@apsara. bare, ernet.in 
(February 1, 2008) 

Abstract 



O 

I The summation of the partial wave series for Coulomb scattering amplitude, 

m . f c (9) is avoided because the series is oscillatorily and divergent. Instead, 

f c {6) is obtained by solving the Schrodinger equation in parabolic cylin- 
>> \ ■ drical co-ordinates which is not a general method. Here, we show that a 

reconstructed series, (1 — cos9) 2 f c (9), is both convergent and analytically 
summable. 
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The partial wave analysis is the most widely used method to obtain cross-sections of the 
scattering processes due to central potentials in atomic, nuclear and molecular physics [1]. 
In this method the scattering matrix Si is the most important theoretical ingredient which 
is the characteristic of a scattering potential. By knowing Si one can obtain the scattering 
amplitude, f{0), and hence the cross-section , a(8)(= \f(6)\ 2 ). The partial wave series for 
the scattering amplitude which sums over discrete angular momentum, I, is written in terms 
of Si as [1] 
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f(0) = —j:(2l + l)[S l -l]P l (cose). (1) 

Despite the importance of an S-matrix it is ironical that the Coulomb potential, V c {r) = 
Q 1 ® 2 } is unique example admitting a simple, exact and analytic S-matrix for all angular 
momenta and energies. The Coulomb S-matrix is expressed as [1] 
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where the super-script C stands for Coulomb, 07 is the Coulomb phase-shift and r/(= Q ^ fc 2M ) 
is the Sommerfeld parameter, T(z) is the gamma function. Surprisingly, in books [1] the 
summation or summability of the series (1) for the Coulomb potential remains not only 
elusive but also without a remark. Generally, the next thing that is done [1] is to take resort 
to solving Schr"0dinger equation in parabolic cylindrical co-ordinates to extract f c (6) as [1] 

f C \0) = - 2ks J {9/2) expHfo log(sin 2 (£/2)) - 2a }] (3) 

In this situation, one may wonder if the partial wave analysis is not a general method 
for obtaining the scattering amplitude for a central potential. Furthermore, if one tries to 
sum the series (1) inserting (2), the exercise turns out to be rather frustrating in that the 
series oscillates showing no convergence, irrespective of number of partial waves included 
[2,3]. Hence, one fails to reproduce f c '(0) (3) even numerically. In fact Eq.(3) serves as the 
benchmark result in the theory of scattering. For the charged particle scattering e.g., proton 
scattering from nucleus where in addition to the Coulomb potential, a short ranged nuclear 
potential, Vjv(r), is used to calculate the scattering amplitude, 

f(°) = ^P 2/ + !)[ ex P( 2 ^) " l]^(cos0), (4) 

where r]i is the phase shift of the combined nuclear and Coulomb potential. Once again, 
the series (4) would show oscillatory divergence and in order to suppress this divergence one 
writes rji — 81 + 07. Notice that 5i is the 07- subtracted phase shift which will also depend on 
the parameter, 77, in some way which, in turn, is not known explicitly. If rji corresponds to 
Vn(t) + Vc(r) and 0/ corresponds to Vcr(r), one can really not assert as to which potential 
form in terms of Vc(r) and Vjy(r), the phase shift 5i would actually and exactly correspond 
to? Only in an approximate calculation such as Born approximation [1] the phase shifts add 
for two potentials. In any case Eq.(4) can be algebraically split as [1] 
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f(°) = E( 2/ + l)[exp(2i<7,) - l]Pi(cos9) + E( 2Z + l)exp(2i<T,)[exp(2i<J,) - l]P,(cos0) 

l%k \l=0 1=0 y 

(5a) 



Subsequently, in the above expression the first series, instead of summing, it is ingeniously 
substituted by f c (9) (3) as 

1 00 

f(6) = f c (9) + — J2(2l + l)exp(2^)[exp(2^) - l]P,(cos0). (5b) 
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By doing so, the convergence of f(9) is solely controlled by the /- dependence of [exp(2iSi) — 1] 
which generally vanishes for large I. In the standard text-books [1] Eq. (5b) is essentially 
written; however, the motivation for doing so is not mentioned. It must be emphasized here 
that the Eq. (5b) does not serve any purpose other than to disentangle the divergence of 
(4, 5a). Equation (5b) is mistaken to separate out the effect of Coulomb interaction from 
f(9), in fact some effect of the Coulomb part goes implicitly in the second part. 

Alternatively, Yennie et al., [4] proposed to bypass this split-up (5b) of f{9) for charged 
particle scattering or for any other instance of scattering where the series for f{9) diverges. 
They suggested an intuitive numerical recipe to calculate f(9) from Eq.(4) itself. According 
to them, one should rather reconstruct a reduced series i.e., (1 — cos9) m f(9) by choosing 
a positive integer, m, such that the resulting series shows convergence. They have success- 
fully calculated cross-sections of high energy electron scattering by nuclei. Somehow, this 
prescription has not received as much attention as it deserves. 

Utilizing this prescription [4] we have been curious to know whether we can overcome the 
divergence problem of Eq.(l) with Eq.(2) and extract the Coulomb scattering amplitude (3) 
directly therefrom. We find that the reconstructed series (1 — cos#) 2 f c (9) is uniformly and 
absolutely convergent and can be summed analytically to recover (3). It may be worthwhile 
to mention that several elegant and rigorous methods of the summation of the partial wave 
amplitudes for long range potentials including the Coulomb potential have been proposed 
[7]. These works are of more general nature which are numerical and which dwell more upon 
approximating the Legendre polynomial Pi(cos9) in various elegant ways. 

The divergence of (1) for Sf can be at once realized by noticing the large / limit of Sf . 
Using an asymptotic property of the gamma function i.e., ll ™ lt F(z + a) — > T(z) exp[a(log z)], 
[5] to find that f™*Sf -> (I + l) 2ir >. It shows that the quantity [Sf - 1] oscillates even for 
asymptotically large values of I. 

Using a property of the Legendre polynomials [5,6] i.e., X^(2Z + l)P/(cos0o)-fH cos 0) = 
2S(cos9 — cos 9), we can write 

f c \9) = (2ik)~ 1 (f^(2l + l)Sf P,(cos0) - 25(1 - costf)^ , (6) 

where S(z) is the Dirac-delta function. This shows an obvious singularity (divergence) in f c 
at 9 = 0. Let us see if we can suppress this divergence. Recalling an interesting property of 
Dirac-delta function i.e., zS(z) = 0, it is tempting to multiply (6) by (1 — cos 9) on both the 
sides. We then get rid of the second term and find 
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(1 - cos6)f c (6) = J2(2l + l)Sf(l - cos^)Pi(cos^). (7) 



Next we reconstruct this series by using a recurrence formula of Legendre polynomials [5,6] 

(21 + l)xP t (x) = (l + l)P l+1 + IP^x), (8) 
where cos 9 = x and we get 

oo 

(1 - x)f(x) = (2ik)~ l + l)S?Pi(x) - (I + l)SfP l+1 (x) - ISfPUx)]. (9) 

o 

oo 

(1 - x)f(x) = (2ik)^ 1 £[(2Z + l)Sf - ZS£x - (I + l)Sf +1 Wx), (10) 

o 

which in terms of gamma functions reads as 

** - -™ - 1 (<» + i)^s - a - * + ^ik^s) *<*>■ 

(ii) 

By multiple usage of the recurrence formula for the gamma function viz., F(z + 1) = ^r(z), 
we can write 

2i*(l - x)f(x) = £ 2^(2/ + 1) ( r^gl^ ) (12) 

For a use in the sequel, here let us denote the quantity appearing inside the large bracket 
in the above equation as 7]. The above series can be re- written in terms of Coulomb phase- 
shifts as 

According to Weierstrass' M-Test [5] , a sufficient condition for a series, {[/ n (2:)}, to be 
uniformly and absolutely convergent is that the series, {M n }, converges. Here, M n > \U n (z)\ 
for the given range of z, where M n must be independent of z. Noting that < 1 for 

— 1 < x < 1, we choose 

= Wj> . (14) 

/ (/ 2 + r ? 2 ){(/ + l) 2 + r ? 2 } 



The series, {M/} upon comparing with X)o° lA (divergent) is divergent. Hence, the M-test 
turns out to be negative for the uniform and absolute convergence of the series (12,13). In 
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fact, the M-test being only sufficient we can not really be assertive about the convergence 
of (12,13). It may be remarked that due to the presence of Pi(x) in Eq.(13) only M-test is 
feasible here. Therefore, we further reduce the series (12) by multiplying by (1 — x) on both 
the sides in anticipation of a better alternative series for which the M-test can be positive 
about its convergence and hence about its summability. The second reduced series with the 
help of Eqs. (8,12) is obtained as 

oo 

2ik{\ - xff{x) = 2 V 2 £[(2Z + mP^x) -(1 + l)T t P l+1 - WP^x)}. (15) 

o 

Carrying out similar manipulations as done earlier from Eqs. (10) to (12), we reconstruct 
the series (15) as 

2i*(l - xff{x) = - i V f £(2/ + l) ^j~3 + ^ P^)- (16) 

Eliminating gamma functions from this series, we get a more transparent expression in terms 
of Coulomb phase-shifts. 

**(! - = -Vd - in? t ^^Xi-^K,-!-^) ^'- < 17 > 

Once again we carry out the M-test, choosing the M-series to be 

M, = t (2 ' + 1> (18) 

o Ml + 2) 2 + + l) 2 + 7? 2 }{Z 2 + 7? 2 }{(/ - l) 2 + 77 2 } 



which upon comparison with X)o° (convergent), is convergent. This establishes the uni- 
form and absolute convergence of (15,16). Let us now sum it up. To this end, we make use 
of a very interesting formula due to Bateman [6] 

(i - *r = z p t 2 l + \t!\ p ^ ^ 

V n + p + 1 (1 + p) n 
where {£) n = r(£ + n)/T{£). By employing Eq.(18) with p — 1 — irj, Eq.(16) yields 

2,:M1 - 4< f( X ) = _^(l-^2-^^)--T(-l + „^ (20) 

By using the identity, zT(z) = T(l + z) and putting x = cos#, we straightaway obtain the 
Coulomb scattering amplitude (3). Thus, the derivation of Eq. (3) presented here supple- 
ments the partial wave analysis for the Coulomb potential. This exercise also demonstrates 
how the divergence of a series is disentangled to extract a physical result. 
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